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(— I I Abstract. Sums of the form X^„<2, E^{n) (fc G N fixed) are investigated, where 



^^(T) = + it)\'' dt - T(iog £ + 27 - 1) 



is the error term in the mean square formula for \C{\ + it)\. The emphasis is on 
the case fc = 1, which is more difficult than the corresponding sum for the divisor 
1^ I problem. The analysis requires bounds for the irrationality measure of e^'^'" and 

^S) ■ for the partial quotients in its continued fraction expansion. 

o 

O ■ 1. Introduction 

Recently J. Furuya [3] investigated the sums of A'^(n) (n G N) when /c = 2 or 
^ ' /c = 3, and as usual 



A(x) = ^ d{n) — x(logx + 27 — 1) 



n<ix 



denotes the error term in the classical Dirichlet divisor problem {d{n) is the number 
of divisors of n, and 7 = — r'(l) = 0.577215649. . . is Euler's constant). For an 
account of the divisor problem, the reader is referred to [7] and [8] . Furuya proved 
that (ci, C2 are suitable constants) 

A'^(n) = / A^(u)dwH — x log^ x + cix log x + C2X + R(x) , 
^ Ji 6 
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where R(x) = 0(x^/^loga;) and at the same time R(x) = Q,±(x^/'^\ogx). Here 
and below, f{x) = Q±{g{x)) means that both 

limsup f{x)/g{x) > 0, \immf f{x)/g{x) < 

hold. This improves much on G.H. Hardy's classical result [5] that R{x) = 
Og(x^"'"'^), where the subscript means that the O-constant depends only on s, an 
arbitrarily small positive constant. The work [3] also brings forth a corresponding 
asymptotic formula for ^„<j. A'^(n), where again the error term is precisely de- 
termined. It seems natural to try to generalize Furuya'a results to sums of £'^(n), 
where A; e N is fixed, ({s) = X^^i '^"^ (3?es > 1) is the Riemann zeta-function, 
and 

(1.1) = / |C(i+^t)pdt-r(log£ + 27-l) 

is the error term in the mean square formula for \({^ + it)\. This is motivated 
by the general analogy between A (a;) and E{T), inspired by the pioneering work 
of F.V. Atkinson [1], who produced an explicit formula for E(T) with a good 
error term (see also [7] and [8]). Later research established more precise analogies 
between the two functions, and in particular we refer to the works of M. Jutila 
[14], [15], and the second author [10], [11]. 

In considering the sums 'Ylin<x-^^^'^)i turns out that already ioi k — 1 the 
problem is somewhat different than the corresponding problem for sums of A(n). 
Namely it reduces to the evaluation of 

(1.2) ^ mm+it)?^t [m=t-[t]-\)- 

To see this, note that 

^ E{n) = J2 (^(^) - tt) + TTX + 0(1) 



n<x n<x 



7rx + 0{l)+ [ {E{u) -TT)d[u] 
Jo 

TTX + 0(1) + G{x) - I {E{u) - tt) dil}{u) 

Jo 

7rx + 0{\ogx) + G{x) + 0{\E{x)\)+ f iP{t)\C{l + it)^ dt, 

Jo 



where 



(1.3) G{x):= r {E{u)-TT)du, G{x) = 0{x^/^), G{x) = n±{x^/^). 
Jo 
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The results on G{x) in (1.3) are to be found in the second author's joint work with 
J.L. Hafner [4]. In what concerns the summatory function of A(n), we recall that 
already G.F. Voronoi [19] proved that 

(1.4) ^A(n) = ia;loga;+(7-i)a; + i?o(a;), 

where Ro{x) = 0{x^/'^) and at the same time Ro{x) = f2±(a;^/^). 
In working out the proof of (1.4), one actually obtains 

(1.5) ^ A(n) = |xlogx + (7- |)a; + A(x) + y A{t) dt + 0(logx), 

n<x ^ 

and (1.4) follows from (1.5) when one recalls Voronoi's classical results [19] that 

(1.6) r A{t)dt=^ + Ri{x), Ri{x) = 0{x^'^), Rt{x) = n±{x^'^). 
Ji 4 

The results in (1.3) on G{x) are in fact the analogues of Voronoi's results in (1.6). 

Let us point out the basic difference between the (discrete) function A(n) and 
the (continuous) function E{n) in the context of their respective summatory func- 
tions. For the former we start (we may clearly assume that a; e N in what follows) 
with 

(1.7) J]A(n)= J]d(m)( ^ l) - J] n(logn + 27 - 1). 

n<x m<x m<n<x n<x 

The last sum is easily evaluated by the Euler-Maclaurin summation formula. Since 
we have 

(1.8) ^ l=a;-m+l (m, X e N), 

m<n<x 

one arrives without any major difficulties at (1.5). But in the analogue of (1.7) 
for E{n), the first sum on the right-hand side of (1.7) will correspond to 

(1.9) ^ICa + ^^)P( E l)^^- 

Observe that in (1.9) t is not an integer (except for 0{x) exceptions), and an 
identity like (1.8) will not hold for the sum over n in (1.9). Instead, the function 
[t] = t — il){t) — I will come into play, producing (1.2). It is rather curious, as 
will be shown in the next section, that our problem will reduce to bounds for the 
irrationality measure of e^'^"^ and for the partial quotients in its continued fraction 
expansion. 

Acknowledgements. We wish to thank Prof. Michel Waldschmidt for valu- 
able remarks. 
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2. The summatory function of E{n) 

In this section we shall continue the discussion on the summatory function of 
E{n). In Section 1 we established that 

(2.1) 5^ E(n) = TTX + / il;{t)\a^+it)\^dt + G{x)+0{x^/^), 

where G{x) is given by (1.3), and where we used the (weak) bound (see [7] or [8]) 
E{x) <^ x^/^. Since ip(t) is an oscillating function, one expects that there will be a 
lot of cancellations in the integral in (2.1), but this turns out to be hard to prove. 

To bound the integral (2.1), we first consider subintegrals over intervals of the 
form [T, 2T] . Then we use the approximate functional equation for C'^{s) (see e.g., 
[7] and [10]) in the form 

(2.2) |C(i+it)r = 2 d{n)n-'/' COS (tlog^-t-^)+n{t), 

n<t/(27r) ^ ^ 

where the error term function 7?.(t) satisfies the mean square relation 



(2.3) 



n^it) dt = ^ Q^y^^b rj.^ (A > 0). 



The asymptotic formula (2.3) was proved by I. Kiuchi [16], who used in his proof 
precise asymptotic expressions for 'R-{t) (related to A(t/(27r)), due to Y. Motohashi 
[17]. Using (2.3) we obtain, by the Cauchy-Schwarz inequality for integrals, 

Jt 

j-2T / ^ ^\ 

+ V d{n)n~^/^ V(^)cos(tlog t--]dt. 

^J^/^ Jmax(T,2,rn) V 4 J 

For the integral on the right-hand side we use the familiar Fourier series 
(2.4) = _i f 

m=l 

which is boundedly convergent and thus may be integrated termwise. This leads 
to 

-2T 



Jt 

2T 

-- ^ — ^ d{n)n~^/'^ I sin(27rmt) cos U log t - dt. 



i<logT n<T/7r 



max(T,27rn) 
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Namely, when we write the sines and cosines as exponentials, the relevant expo- 
nential integrals will be of the form 

2T ^ T 



e±^^(*)dt l<n<- 

max(T,27rn) V 



where 



F(t) = tlog t - 27rmt, F'(t) = log 27rm, F"(t) = -. 

zTrn zTrn t 

Therefore for m > logT we have |-F'(t)| ^ m, hence by the first derivative test 
for exponential integrals (see e.g., Lemma 2.1 of [7]) the contribution of such m is 
0(TV2iogT). Next, for T/{2n) <n< T/tt, we have 

\F'(t)\ = 27rm - log — ^ > 27rm - log > 27rm - log2 > 7rm (m e N), 

hence the contribution of such n is again 0{T^^'^logT). Likewise the contribution 
of satisfying m > (l/(27r) + s) logT is 0(T3/4). ^j^^^^ 

remains is 

(2.5) 

f.2T 

V'(^)|C(i + ^^)rd^ = OiT^/') 

2T 

-- V — V / sin(27rmt) COS I Uog ^ t - - ) dt. 

m<(l/(27r)+e)logT n<T/(27r) T 

The saddle point (the root of F'{t) = 0) is to = 27rne2'^"^, and it lies in [T, 2T] 
for n X Te~^'^"^. Using the saddle-point method (direct integration over a suitable 
contour in the complex plane joining the points T, 2T and passing through to) we 
obtain that the expression in (2.5) is asymptotic to sums of the type (C > is a 
generic constant) 

(2.6) C — Yl exp (-27rme2'^"^) . 

m<(l/(27r)-|-e)logT nxTe-^'^'" 

The problem is thus reduced to the estimation of the sum over n in (2.6). One 
can use the results of a classic paper by J.R. Wilton [21]. He proved, among other 
things, the functional equation 

(2.7) D{x, rj) = T]-^D{7]^x, -rj'^) + 0{x^/^ logx), D{x, v) = d{n)e^''''^'' , 

n<x 
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where rf'x » 1 and < 77 < 1 is real. After taking the conjugate sum in (2.6), (2.7) 
can be applied with x = CTe"^'^'^, r] = {e^'^"^}, where {y} = y — [y] denotes the 
fractional part of y. The problem is that it is difficult to find an explicit expression 
for T] — r]{m). On the other hand, Wilton (op. cit.) states that for irrational rj 
one always has 

(2.8) D{x,r]) = o(a;loga;) {x ^ 00) 

uniformly mrj — r]{m). The bound (2.8) can be applied in our case, since the num- 
bers r]{m) are all transcendental. Namely, since the number e'^ is transcendental, 
all the numbers e^"^"^ are also transcendental. The fact that — (— 1)~* is tran- 
scendental follows from the solution of Hilbert's 7th problem (see e.g., C.L. Siegel 
[18] for a proof) that is transcendental for algebraic a 7^ 0, 1 and irrational (3. 
This was proved by Gel'fond and Schneider independently in 1934. Wilton [21] 
also proved the following: let 

1 11 

- = ao + — ■ ■ • • < 7/ < 1 

rj ai+a2+ 

be the expansion of l/?7 as a simple continued fraction. If rj is irrational and the 
convergents an satisfy <C n^+^ {K >0), then 



2+K 

X, 



(2.9) ^ d{n) sin(27rn?7) < x^/^ log 

n<x 

while if 

(2.10) an < e^^ {K > 0) 
holds, then 

(2.11) Yl exp(27rmr7) « logx (h = . 

Therefore the information on the size of the a^'s would lead, by (2.9) or (2.11), 
to good bounds for ^„<^ d{n) sin(27m?7), which is the sum we need to estimate. 
Additional care should be exerted because in our case r] = e^'^'", so that = 
an{iTi) will depend on m as well. To determine the true order of = cin{'m) 
seems to be a difficult problem; even the fact that e'^ is not a Liouville number is 
not known. In the next section we shall provide an explicit bound for the partial 
quotients of c^^"^. Although poorer than (2.10), it is still sufficient to provide a 
non-trivial bound which in our case improves on the general bound (2.8) in the 
case when rj is irrational. 
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3. The results on the sum of E{n) 

We begin by recalling a bound for the irrationality measure of e'^"^. Namely, 
for any positive integers m, p and q with p > 3, we have 



(3.1) 



> 



exp I — 2^^ (log 2m) (log p) (log log p) | . 



This result was proved by M. Waldschmidt [20, p. 473]. Now for any non-zero 
integer m let 

e'^^ = [ao(m);oi(m),...] 
be the expansion of e'^"^ as a continued fraction. 

LEMMA 1. There exists an absolute positive constant c such that 
(3.2) log log Gn (m) < c(n + log m) log(n + log m) 

holds for any non-zero integer m and any non-negative integer n. 

It is sufficient to establish Lemma 1 for positive values of m, since e~'^'^ = 
[0; ao(m),ai(m), . . .] for m > 1. Let ^ be a positive real number with continued 
fraction expansion 

^ = [ao;ai,a2,...] 

and set, for n > 1, 

— = [ao; ai, a2, . . . , On]. 

From the theory of continued fractions it is known that 

1 



(3.3) 

It follows from (3.1) that 



C - — 



< 



(3.4) 



QnQn+1 



> exp{ —2^^ (log 2m) (log q) (log log q) } , 



(n>l). 



for all p/q with q > e'^"^, since p < 2qe'^^, giving 

logp < log 2 + log q + 7rm < log 2 + 2 log g < 3 log q. 

In particular, (3.4) is true for the convergents ph{m)/qh{m) to e'^'^ as soon as 
(Ihijn) > e^"^. Let ?i > 2 be an integer. Writing simply Pn/qn instead of 
Pnijn) / InipT) 1 it then follows from (3.3) and (3.4) that 



(3.5) 



qn < exp|2'^^(log2m)(logg^_i)(loglogQ^_i)|, 
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if Qn-i > e^"^. We infer from (3.5) that there exists an absolute positive constant 
ci such that 

(3.6) log2 Qn <Cl+ log2 m + log2 Qn-l + logg Qn-l, 

where, for j > 2, log^- x = log(logj_i x) (log^ x = logx) denotes the j-th. iteration 
of the logarithm. On dividing (3.6) by log3 qn, we obtain 

(3.7) < 1 + (ci + log2m) + -j . 

logs Qn logs gn logg Qn 

Using repeatedly the fact that Qn > Qn-i, it follows that 

< n + (ci + log2 m) + 



^3 logs Qn logs Qj log3 Qh+l 

ci + log2 m 

< n + n — hlogag/i, 

logs Qh 

where h is the integer such that qh-i < e^'^ and qh > e^'^. This choice of h and 
(3.1) imply that 

(3.9) Qh < exp{c2m(log2m)^}, 

where C2 (as C3, C4 and C5 below) is an absolute positive constant. Consequently, 
we infer from (3.8) and (3.9) that 

logo Qn Ci + logo m , / 1 \ 

< n + n — h C3 + logm < cAn + logm). 

loga qn log2 TO 

This gives at once the asserted upper bound (3.2), namely 

log2 an < log2 qn < C5(n + logm) log(n + logm). 
Now we rewrite (3.2) as 

(3.10) anim) < exp |(n + logm)'=(^+i°sm) | 

and proceed to estimate, in view of the transformation formula (2.7), 
(3.11) 

D{x,'r]{m)) := ^ (i(n) exp(27ri?7(m)) \r]{m) = e"^"^"*, 1 < m < — - — j. 

n<x ^ 
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Then l/r]{m) = e^"^"*, and the above discussion apphes with m replaced by 2m. 
The above range for m suffices, since by trivial estimation the contribution of 
m > loga T/ logg T to (2.6) is 

•C > • logT <^ riogrexp< — 6 —>. 

m e^'^^ ^ 6 F| logo T J 

m>log2T/log3T 

From Wilton's bounds [21, eqs. (6.3///) and (6.3/v)] we have 

(3.12) D{x,r]{m)) <^ x^^^ log^ a; + min^OAr(2m)a;^/^ logx, 2~^^a;loga;^, 

for large N satisfying <^ log a;. We use (3.10) to deduce that, for 



N < 



log2X 



. c / loga X 

we have 

N 1/9, 1/9, f/",. logoXXcC-'V+logax/logga;)^ 

ajv(2m)a;^/^loga; < x^/^ log a; exp < iV + > 

tV log^xJ J 

<a;loga;exp f-c|^^^ (C = C(c) > 0), 

while for A'' > (l/c + o(l)) log2 x/ log^ x obviously 

2~2^a;loga; <C xlogxexp ( —C- — - — 

V logs X 

This means that we obtain the following non-trivial result, which improves on 
(2.8), although only by a quantity that is less than (logx)^ for any given e > 0. 



THEOREM 1. If D{x,r]{m)) is given by (3.11), then for 1 < m < and 
some constant C > we have 

(3.13) D(x^r]{m)) -C a; log a; exp (—C-—^^ — 

V logs X 



To apply Theorem 1 we transform the sum over n in (2.6) by (2.7) (in the 
range 1 < m < logaT/loggT) with x = CTe-^'^"^, r/ = e^^"". We obtain 
D{x^ rj{m)) with rj{m) = e~^'^"^, x = CTe^'^"^, plus an error term of 0(T^/^ logT). 
We use (3.13) for this D{x,r]{m)), estimate the contribution of the range m > 
log2 T/ log3 T trivially, and arrive at 
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(3.15) H{x) <^ xlogxexpl— C 



THEOREM 2. We have 
(3.14) XI ^(^) = TTX + H(x), 

n<x 

where, for some C > 0, unconditionally 

If (2.9) holds with rj = e^'^^ uniformly for 1 < m <^ log a;, then we have 
(3.16) H{x) = Oix^/"^). 

Remark 1. The intrinsic difficulty in this problem is to bound the sum 
D{x,rj{m)). In fact, unconditionally our bound for H{x) in (3.15) is poorer than 
the main term ttx in (3.14). However, if the bound in (2.9) is true in our case, 
then (3.16) can be complemented with H{x) — fl±{x^/'^). Namely instead of using 
(2.2)-(2.3), we may use an (unsymmetric) approximate functional equation for 
+ it)\'^ with a sharp error term (see e.g., [7, Chapter 4] ). In fact it may be 
conjectured that instead of (3.15) with (3.16) one has 

E{n) = 7rx + G{x) + 0e{x^/^+'), 



n<x 



where G{x) is given by (1.3). 

Remark 2. Another discrete sum involving the divisor function was investi- 
gated by Coppola-Salerno [2] and A. Ivic [12]. The former have shown that, for 
<h< |v^, L = logT, 

(3.17) (^(^ + ^) - ^(^^))' = ^TUXog"" (^) + 0{TUL'>'''Vl). 

T<n<2T ^ 

This is an asymptotic formula with a weak error term, improved in [12] by the 
second author to 

E (A(n + C/)-A(n))^=TC/^c,log^(^) 

(O.lo) T<n<2T 3=0 

where 1 -C C/ < ^VT with cs = Stt"^. However, the analogue of (3.17) or (3.18) 
for the sum 

(3.19) Yl {E{n + U)-E{n)f 

T<n<2T 

does not carry over, for the same reason for which we had difficulties in evaluating 
sums of E{n)-. because E{T) is a continuous function, while A(a:) is not, having 
jumps at natural numbers of order at most Os{x^)- The true order of magnitude 
of the sum in (3.19) seems elusive. 
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4. Sums of E''{n) for A; > 1 
We begin by noting that in general, for k fixed, one has 



(4.1) Y^E^{n)= f E\u)d[u]^ [ E^{u)du- f E^ 

-^0 Jo Jo 



(u) d'i/j{u). 



n<x 

Integration by parts gives 



(4.2) 



I E''{u)d^{u) = E''{x)^{x)-k ij{u)E''-\u)E'{u)du 
Jo 

« 1 + \E{x)\'^ + £ \E{u)\''-' (|C(i + iu)\' + log(|-) + 27) du, 

where we used the defining property (1.1). Here we distinguish two cases. 

If k = 2, 4, then by the Cauchy-Schwarz inequahty for integrals the last integral 
above is 

« E^''-\u)du- £(\C{l+iu)\''+log''u^du^ ' 
« x'^/^ + [x^^h^'^-^) . X log^ x) « x""!^ + log2 ^ 

where we used (see the works of D.R. Heath-Brown [6], Ivic-Sargos [13] and W. 
Zhai [22]), 



(4.3) 



/ \E{u)\^du ^ x^"^^l^ (A; e N, 1 < /c < 9). 
Jo 



Larger even values of k could be, of course, also considered, but for such values 
we do not have precise formulas for the integrals of as we have when A; = 2 or 
A; = 4. Namely we have (see [8]) 

(4.4) / E'^ {u)du = C^x^l'^ + O (a; log^ x) (C2 > 0) , 
Jo 

and (see W. Zhai [22, Part III]) 

(4.5) E\u) du = C4x'^ + Os(a;2-3/28+e) (c^ > q). 
Jo 
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If k = 3,5,7,9 then A; - 1 is even, and we have \E{u)\''~'^ = E''~''-{u). Since 
(4.3) holds and E{x) <S x^/^, we see that the last integral in (4.2) is 

/ E''-\u)(e'{u) + 0{logu)) du <^ \E{x)\'' + x^''+^^/^logx <^ x''/^logx. 
Jio ^ 

Note that we also have (see W. Zhai [23], who improved the exponent 7/4—1/12 
of Ivic-Sargos [13] in the error term) 

(4.6) ^ 

E^u) du = CW^ + 0,(a;^/^-S3/393+s) ^ g ^lligs . . . , C3 > , 

and, in general (see W. Zhai [22]) for 5 < A; < 9, 

(4.7) r E^{u) du = Ckx^+^'^ + 0{x^^^/^-^^^^) {Ck > 0), 
Jo 

where S{k) is a positive constant which may be explicitly evaluated. From (4.1)- 
(4.7) it follows that we have proved 

THEOREM 3. We have 



J2 E\n) = C^x^l'' + 0(a;^/4log'a;), 

n<x 

^ E\n) = Csx'/' + Os{xy^-''/^^^+'), 

n<x 

Y,E\n) = C,x' + 0,{x'-'/''+% 

n<x 

and for k = 5,7, 9, 

J2 E\n) = + C>(a;^+'=/^-'^(^)) + 0{x''/^ logx). 

n<x 

Naturally, other values of k could be also considered, but as already remarked 
the results would not be very good in view of the existing results on E^{u) du. 
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